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A massive exoplanet candidate around KOI- 13: Independent 
confirmation by ellipsoidal variations 
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ABSTRACT 

We present an analysis of the KOI- 13. 01 candidate exoplanet system included in the Septem- 
ber 201 1 Kepler data release. The host star is a known and relatively bright (otkp = 9.95) vi- 
sual binary with a separation significantly smaller (0.8 arcsec) than the size of a Kepler pixel (4 
arcsec per pixel). The Kepler light curve shows both primary and secondary eclipses, as well 
as significant out-of-eclipse light curve variations. We confirm that the transit occurs round 
the brighter of the two stars. We model the relative contributions from (i) thermal emission 
from the companion, (ii) planetary reflected light, (iii) doppler beaming, and (iv) ellipsoidal 
variations in the host-star arising from the tidal distortion of the host star by its companion. 
Our analysis, based on the light curve alone, enables us to constrain the mass of the KOI- 13. 01 
companion to be Mc - 8.3 + 1.25Mj and thus demonstrates that the transiting companion is 
a planet (rather than a brown dwarf which was recently proposed by Szabo (201 1)). The high 
temperature of the host star (Spectral Type A5-7V, T^s - 8511 - 8020 K), combined with 
the proximity of its companion KOI- 13.01, may make it one of the hottest exoplanets known, 
with a detectable thermal contribution to the light curve even in the Kepler optical passband. 
However, the single passband of the Kepler light curve does not enable us to unambiguously 
distinguish between the thermal and reflected components of the planetary emission. Infrared 
observations may help to break the degeneracy, while radial velocity follow-up with cr ~ 100 
m s ' precision should confirm the mass of the planet. 
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1 INTRODUCTION 

In Mislis et al. (2011) we demonstrated the potential for estimating 
planetary masses solely from light curve (LC) data for a restricted 
sample of systems. The distinguishing characteristic of these sys- 
tems is that the planet must be both massive enough and close 
enough to the host star to induce significant tidal distortions in the 
stellar ellipsoid. Such non-sphericity in the plane of the sky for 
a rotating star will lead to a periodic photometric signal which is 
detectable given sufficient signal-to-noise. In Mislis et al. (2011) 
we applied the technique to the transiting exoplanet HAT-P-7 ob- 
served with Kepler, and estimated the mass of the planet to be 
Mp(LC) = 1.27 ± 0.12My which is very close to the RV solution, 
Mp(RV) = 1.20 ± 0.05 (Welsh et al. 2010). 

The second exoplanet candidate list released by t he Kepler 
team in February 201 1, and published in l iBoruckillioilh . contains 
1235 exoplanet candidates. The light curves span measurements 
taken between May 2009 and September 2011. Motivated by our 
previous work, we searched the light curves of all 1235 planetary 
candidates for evidence of ellipsoidal variation. The brightest and 
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best candidate for a star showing strong out-of-eclipse ellipsoidal 
variations is KOI-13.0]Q candidate, which is identified with an A5- 
7V type star (with a close companion 0.8 arcseconds to the wesQ). 
The light curve magnitude is otkp = 9.95 (Kepler passband), and 
shows primary and secondary eclipses with a published period of 
1.7635892 days. The shape of the eclipses are indicative of either a 
small transiting companion, or a blended system, whereby a deeper 
eclipse is diluted by contaminating light from nearby stars in the 
same Kepler pixels. 

KOI-13.01 is already the subject of a through analysis of the 
source characteristics and the transit shape by Szabo et al. (2011). 
They conclude that the system comprises a double A-star wide 
binary system, the more massive component of which is being 
eclipsed by a brown-dwarf or very-low-mass star on the basis of 
a radius determination = 2.2 ± 0.3Rj (Rc' companion radius). 
In this paper we independently confirm that the brighter compo- 
nent of the visual binary is the host star for the eclipses (Section|2]l. 

' Shortly before submission of this manuscript we became aware of an 
analysis of the same system, drawing similar conclusions, submitted by 
Mazeh et al. (201 larXivl 1 10.3512). 
2 Szabo etal. 2011 
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Barnes et al. (2011) use a more sophisticated analysis to measure 
the radius of the host star and transiting companion, including the 
effects of a gravity-darkened rapidly-rotating host star. The updated 
values (Rc = 1.4 ± 0m6R,,R^, = 1.77 ± O.OURe) are significantly 
smaller than those found by Szabo et al. 201 1 , and lead Barnes et al. 
(201 1) to conclude that KOI- 13.01 may be a planet. Finally, Rowe 
et al. (2011) presented the system at the 2011 AAS meeting and 
also suggested the companion is a planet. 

In Section |3] we perform modeling of the out-of-eclipse pho- 
tometric variation, which puts very strong constraints on the mass 
of the eclipsing/transiting companion, and lends support to the ar- 
gument that KOI13.01 has a planetary mass. 



2 KEPLER DATA 

We have used the Kepler public light curve and centroid curve 
data available from the Kepler archive hosted by the Multimission 
Archive at STScI (MASTUfor our analysis. The observations com- 
prised four long-cadence and four short-cadence datasets (spanning 
Kepler Q1-Q3). The data were processed by the Pre-Search Data 
Conditioning (PDC) module (Kepler Data Processing Handbookfl 
For objects brighter than otkp = 11.3 some saturation, and then 
bleeding, of the CCD pixels will occur. This charge is not lost, 
and using an appropriately shaped mask will preserve the flux, 
and maintain precise photometry even for extremely bright sources 
(Gilliland et al. 2010, 2011). Both short-cadence and long-cadence 
data are equally affected given that they use the same 6 second ex- 
posure. Figure[T]we show the Kepler light curve folded on the pe- 
riod of 1.7635892 days found by Borucki et al. (201 1), with phase 
zero set to be the center of the primary minimum. 

2.1 Identification of the companion host star 

Borucki et al. (2011) suggest that the source associated with KOI- 
13 is actually a double star (separation 0.8 arcsec, Amag = 0.4). 
Szabo et al. (2011) show that the Kepler source is coincident with 
the double star (BD+46 2629A and B, hereafter K0I-13A & KOI- 
13B), two stars of spectral class A, separated by 1.18±0.02 arcsec- 
onds with a magnitude difference of Amag = 0.20 ± 0.04 magni- 
tudes. 

Howell et al. (2011) use speckle imaging to find a separation 
of 1.16 arcseconds and Amag = 0.85 at 562nm (Amag = 1.03 at 
692nm) between the two components. However, this photometric 
difference between the two stars is in serious disagreement with 
the Szabo et al. (2011) photometry. Having inspected the Szabo et 
al. (2011) image, we have opted to adopt their value for the mag- 
nitude difference between the two stars in the Kepler bandpass (as 
do Barnes et al. 201 1). We note that this is a V-band measurement, 
which is not identical to the Kepler passband, but for an A-type 
star, the differences are small. In our analysis (Section 5) we con- 
sider the effects of changing the contribution of this companion star 
to the Kepler light curve. Based on an unresolved spectrum, and a 
combination of resolved and unresolved photometry, Szabo et al. 
use model isochrones to find effective temperatures for the A and 
B components of T^^ = 8511K (8128~-8913K) and Tf^ = 8222K 
(7852~8610K) respectively. 

The centroid curves for KOI- 13 are shown in the middle panel 

' http://archive.stsci.edu 
available via |http://keplergo. arc.nasa.gov/Documentation. shtml| 
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Figure 1. Phase folded light curve and the secondary eclipse of KOI-13 
system. The primary transit has been removed to enhance the visibility of 
the out-of-transit variabiUty. 

of Figure[2] Batalha et al. (2010) have demonstrated that the very 
tiny shifts in centroid, expected for blended occulting sources, can 
be measured accurately from the Kepler data. Thus they can be 
used to show where the largest variations in flux are occurring with 
respect to the center-of-light. The strong periodic offset in the Ke- 
pler centroid curve data agrees extremely well with the time of pri- 
mary transit. Figure [2] (bottom) also shows an illustration of the 
direction (and magnified amplitude) of the source centroid during 
primary eclipse. This can be explained simply, if K0I-13A be- 
comes fainter during every transit and the center-of-light moves 
towards the west (towards KOI-13B). This interpretation indepen- 
dently confirms the conclusions from Szabo et al. (201 1) who used 
both the Kepler pixel light curves, and additional ground-based 
lucky imaging. Thus we confirm that KOI-13A is the host to the 
companion. 

2.2 Removal of the blended light 

For our analysis, the Ql, Q2 and Q3 datasets were used, corre- 
sponding to a light curve spanning 418 days in total. The light curve 
was first corrected to remove the flux from the secondary which we 
assume to be constant. A magnitude difference of Amag = 0.2 was 
used, and 45 % of the total flux is subtracted from the Kepler light 
curve ( iSzabdl201lh . 

The light cur ve was phase-f olded on the best-fit Kepler period 
of P= 1.7635892 teoruckill20Tlh . rebinned by a factor 1770, and 
phase=0 set to be the midpoint of the primary transit. We use the 
updated values for Rc, Rt,, a (semi-major axis of the orbit) and 
inclination (i) from Barnes et al. 2011 (R^, = 1.76 ± 0.014 i^o, /^c = 
1.4 ± 0.014 Rj, a = 0.0367 AU & i=85.9). The phase-folded light 
curve, including the secondary eclipse (but excluding the primary 
transit to emphasize the out-of-eclipse variation), is shown in Fig.[T] 
Finally we found no evidence for additional periods in the Kepler 
corrected light curve, such as Mazeh et al. 2011 suggested (P = 
1.0595 days). 



3 MODELING 

In our modeling, we consider four main components which con- 
tribute to the out-of-eclipse phase dependence of the light curve of 
KOI-13, summarised in Equation[T] 



A massive exoplanet candidate around KOI-13 3 




0.8 
0.6 

S 0.4 

u 

o 

^ 0-2 

OJ 

If) 

B 0.0 

o 

LiJ 

Q -0.2 
-0.4 



K0I-13B 



P. A = 280.0 deq. 



-0.8 -0.6 -0.4 -0.2 0.0 
RA Offset [arcsec] 



0.2 



Figure 2. Top: The transit phase folded hght curve. Middle: X & Y axis 
offset of the centroid. The offset is periodic (same period as the transit) in 
the same phase. Offset variations show that the primary star becomes fainter 
during each transit event. Bottom: RA & Dec coordinates of the KOI-13A 
(right) and KOI-13B (left) plus the RA & Dec offset vectors. The (0,0) point 
refers to the K01-13A coordinates (RA = 286.971° & Dec = 46.8684°). 
The vectors have been magnified by a factor of 100. 



Figure 3. Top: Phase folded light curve and the best fit model (model - 3) 
including all three coinponents (/((,, f[,.{, fc\\,fdop)- The blue solid line refers 
to the /th+/ref component. The additional flux froin ellipsoidal variations at 
phase a = 0.5 is zero, that's why the total flux at the same phase becomes 
only from the /th+/ref component (ignoring the secondary eclipse). Bot- 
tom: The ellipsoidal variation residuals and the best fit (red solid line). The 
inass of the body which causes these variations must be Mq ~ 8.3 Mj. The 
blue solid line refers to a brown dwarf (BD) ellipsoidal model with mass 
Mbd = 13.0Mj. 



stellar spectral energy distribution with respect to the bandpass of 
the instrument (Loeb & Gaudi 2003, Mazeh & Faigler 2010). This 
final component contains no additional parameters. The first four 
components are shown in equations 2-5 (for their derivation, and a 
detailed discussion, see Mislis et al. (2011)). In Mislis et al. (2011), 
we show that if we can disentangle the ellipsoidal variations from 
the close companion's thermal and reflected light components, we 
can solve for the mass of the companion. 



(2) 



ftot 

T 



= 1 + 



fthJay + fth ,night + fref + fell + fdop 



(1) 



In our simple model we include (1) thermal emission from the com- 
panion (day and night side - f,hjiay, fth,night )> (2) reflected light from 
the surface of the companion, (3) ellipsoidal variations due to tidal 
forces between the unseeen companion and the star, and (4) flux 
variations arising from the Doppler shift (Doppler boosting) of the 



- ° ' ' sm (;) 



— — = (3 -p)— 



(3) 



(4) 



(5) 



4 D. Mislis 



where f±l f* is the relative flux from the thermal emission (for the 
day side, /th.day//*, A{e) = (5(6) and for the night side, /,h,n,ght//*, 
A(d) = 1 - O(0)), /rcf//* is the relative flux from the reflected light, 
/ell//* is the relative flux from the ellipsoidal variations and /*,,,//* 
is the Doppler Boosting relative flux. O(0) is the phase function, a 
is the semi-major axis, Const{Tcff,A,a^„\,£) is a constant function 
of the effective temperature of the star and the temperature of the 
planet (Mislis et al. 2011), and [} is the gravity darkening. Finally, 
K is the radial velocity amplitude, p is a function of p{A, T^ff). 

ghc!kAT,ff^2 - hc/kAT,ff) - 3 
P = e'ciUT.ff - 1 

Our model for the flux changes arising from the ellipsoidal distor- 
tion (Eq.O is a simple approximation, and assumes that the system 
is in tidal equilibrium (Rowe et al. 2010, Carter et al. 2011, Mis- 
lis et al. 2011). The aboi and e parameters refer to the bolometric 
albedo and the energy circulation respectively. Both equations[2]& 
[3]are a function of (t>(d), but the mass of the secondary body is a de- 
pendent parameter only in eq.|4]&[5] We note that in our equations 
the terms R-^/a and Rc/cr both appear. These are measured directly 
from the depth and duration of the primary transit and the period of 
the system. For consistancy these values are taken from Barnes et 
al. 2011. 

If we were to solve for and remove the first two components 
(thermal emission & reflected light), the residual light curve will 
contains the mass information (Fig. [3]- bottom). An important re- 
sult of the model is that the ellipsoidal component contributes zero 
to the total light curve flux (assuming e=0) at phase=0.5 (centre of 
the secondary eclipse), and that the phase function depends only 
very slightly on the balance between the day- and night-side prop- 
erties of any companion's emission (i.e. it is extremely difficult to 
measure this with single passband optical data). Thus the thermal 
emission and the reflected light components will be very difficult to 
distinguish in a single filter. Nevertheless, one of the goals of mod- 
eling the full Kepler light curve for KOI- 13 is to try to constrain the 
brightness temperature of KOI-13.01. The main different between 
our model and other algorithms (such as BEER - Faigler & Mazeh 
201 1) are that (i) we investigate two phase functions (Lambert and 
geometrical spheres - Mislis et al. 201 1) and (ii) we include distinct 
reflected and thermal components (including day & night side, Oboi, 
6). The choice of phase functions can alter the mass of the planet 
by 36% (see Section 5). 

KOI-13A is so hot {T^ff = 8511 + 400 K), and the candi- 
date planet so close (a = 0.0367 AU, ~ 4.5/?*), that the thermal 
flux of the candidate is likely to be non-negligible, even at opti- 
cal wavelengths. The equilibrium temperature for any zero-albedo 
companion to KOI-13.01 is T^^ = T.ff^R.jla = 2864K (for 
(R, = 1.694 ± 0.013 Rq, Barnes et al. 201 1), assuming that it does 
not have its own internal source of heating (as one would expect 
for a brown dwarf). The hottest known planet discovered to date is 
probably WASP-33b, with an equilibrium tempera ture of 2750K , 
and a brightness temperature of Tb = "i^lQ^^K ( SmithI ( 1201 ih ). 
measured from ground-based infrared {A = 0.9fim) observations of 
the secondary eclipse (depth of 0.109 per cent). The relative ther- 
mal emission for WASP-33b in the optical (Kepler passband) would 
be 10"* X / (assuming a blackbody). We note that the secondary 
eclipse of KOI-13.01 is 1.2 • lO"" x / (Szabo et al. 2011) in the 
Kepler passband. 

As a starting point, we consider three simple scenarios for the 
phase light curve data alone, i.e. excluding all datapoints within the 
primary transit or the secondary eclipse. In case C,i,enn we assume 
that the bolometric albedo of the companion is fixed at Oboi = 1 -0, 



which means that the planet is perfectly reflective, that the ther- 
mal emission of the planet is zero, and therefore the model con- 
tains only reflected light and ellipsoidal variations. In case C„/;, 
we assume that the geometric albedo is zero (a^ = 0.0), so the 
reflected light is zero and the model contains only thermal emis- 
sion ellipsoidal variations, and Doppler boosting (/j^p is a function 
of p(A, Tgff) and K - Eq. Finally, in case Ci,ybr we consider a 
hybrid case, and allow cig and Oboi to be free parameters, i.e. both 
thermal emission and reflected light are present in the light curve in 
addition to any ellipsoidal variations. So in C„// we fit for tjg and 
Mc, in C therm We fit flboi, e (energy circulation) and Mc, and in C;„6r 
we fit all the parameters above. Eccentricity (e) and the periastron 
angle (w) are free parameters in all three scenarios. 

As a final step, we add the primary transit and secondary 
eclipse data back into the light-curve, and re-assess C;,,*, (and name 
it C*^,;,^), primarily to see if the additional information on the depth 
of the secondary eclipse can provide useful constraints on the ther- 
mal emission form the planet (as has been suggested by Mazeh et 
al. 2011). 

An important and unresolved factor is how we treat the phase 
function of the reflected light. We consider two approaches: (1) 
modelling the planet as a Lambert sphere (Russell 1916), assum- 
ing that the intensity of the reflected light is /„/ = 1 /;r at at phase 
z = n/2 (as in the case of Venus); (2) An alternative choice would 
be to assume that the reflected light is f„f = 0.5 at phase z = n/2. 
In Mislis et al. (2011) we demonstrate that this choice leads to a 
significantly different phase-contribution for the reflected light, and 
therefore has a significant impact on the derived mass of the com- 
panion. The companion mass is significantly larger (by 36 per cent) 
if we use the Lambert sphere compared to the geometrical sphere. 
In our analysis we are unable to distinguish between the two cases. 
Thus, in the next section we present the results from the Lambert 
sphere case, which gives the larger companion mass. 



4 RESULTS 

With all three cases we find that the orbit is circular and that the 
eccentricity value is e = 0.0±0.05. All three scenarios require 
strong ellipsoidal variations to explain the phase light curve. The 
hybrid model Ci,ybr is significantly preferred with a confidence of 
99 per cent confidence (2.55cr) compared to C,i,„m, and 99.6 per 
cent (2.85cr) compared to C„//. Scenarios Ctherm and Crefi are in- 
distinguishable at the Icr level. Table|4]shows the results of all three 
cases (excluding primary transit and secondary eclipse) and Fig.|3] 
(top) shows the Kepler light curve and the components for the Ci,yhr 
scenario. In all three cases, we agree on the mass of the unseen 
companion to within a few per cent, suggesting that the technique is 
largely independent of the precise nature (be it thermal or reflected) 
of the emission arising from the companion. Although note that 
the mass constraint improves as the model becomes more complex. 
Both Ctherm and Crefi are significantly less good than the hybrid 
model, suggesting that components of both thermal and reflected 
light are contributing to the phase light curve. Note that we do not 
give errors on e, Oboi and in Chybr because they are essentially 
unconstrained at the I-cr level, with large degeneracies between flg, 
flboi, and e. 

The mass derived from the best fitting model of the candidate 
is Mc = 8.3 ± 0.9Mj, rather more massive than a typical hot Jupiter, 
but significantly below the brown dwarf or M dwarf mass proposed 
by Szabo et al. (2011). This error on Mc is the error from fitting 
the model, and does not allow for the uncertainty in the mass of 
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Crefl 


Ciherm 


Chybr 




[Reflected] 


[Thermal] 


[Mixed] 


Mc [Mj] 


8.0±1.1 


8.3±1.0 


8.3±0.9 


e 




0.17±0.17 


0.10 


Obol 




0.01+0.1 


0.80 


«g 


0.28±0.1 




0.28±0.28 


e [deg] 


0.0±0.05 


0.0±0.05 


0.0±0.05 


xl 


1.042 


1.032 


1.03 



Table 1. Model results for the three main scenarios discussed in the text. 
From left-to-right the models increase in complexity. We hst the fitted pa- 
rameters and their 1-cr errors. In each case we are fitting 90 datapoints in the 
phase-folded light curve. For scenario Chxbr, the parameters shown without 
errors are essentially unconstrained by the data. 



This then is our motivation for scenario C? , introduced in the 

hybr 

previous Section. We added back into the light curve an additional 
36 (in eclipse) data points, while adding no additional parameters to 
the model, i.e. the shape of the secondary eclipse is completely de- 
scribed by Chybr- A little surprisingly, the additional data do nothing 
to improve our constraints on the nature of the secondary as param- 
eterized by flg, fflboi and e. Perhaps this should not be so surprising, 
the depth of the secondary eclipse in a single filter is unlikely to 
tell us very much about the ratio of thermal to reflected emission 
from the planets atmosphere. Infrared measurements could perhaps 
make a significant impact here, as they could be combined with the 
optical passband to shed light on the colour of the planet. We also 
note that our solution to the mass of KOI- 13.01 is not affected by 
the additional data points (see Figure|4ll. 
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Figure 4. x\ vs Mc for model fits with and without the secondary eclipse 
(Chybr, Cjjjj,^) (solid and dashed curves respectively. The horizontal line il- 
lustrates the l-o" error level for each model. 

the primary, rather this should be seen as a 20 per cent error on the 
mass ratio, Mc/M,=0.0039. We take the mass from Szabo et al. 
(2011) to be 2.05Mo and assume an error of 10 per cent. We also 
consider a 0. 1 per cent error in the radius ratio RqIR, (Barnes et al. 
2011). Accounting for both of these leads to a final constraint on 
the mass of KOI-13.01, Mq = 8.3 ± 1.25Mj. 

4.1 How hot is KOI-13.01? 

Any attempt to constrain the brightness temperature of the compan- 
ion is severely hampered by the problem of distinguishing between 
the reflected and thermal components of our model, exacerbated by 
the single passband. In principle, if there were no reflected light, 
nor ellipsoidal variations, we could use the difference between the 
depth of the secondary eclipse (star only) and the point just before 
ingress to (or after egress from) the primary transit (i.e. star plus 
planet night-side) for an orbit with zero eccentricity (e.g. TrES-2, 
[o'Donovan (2010)). In practice, this is rather more complicated, 
because very small contributions from the ellipsoidal variations, 
and the marginally visible reflected light components are contribut- 
ing to the total light of the system at this point in phase. And they 
need to be included in the modelling to constrain the night-side 
contribution to the system flux. Mazeh et al. (2011) have ignored 
these effects, and estimate the temperature of the companion to 
be 2600K. Our rather more complete treatment of the light curve 
includes both thermal and reflected components, and includes the 
secondary eclipse in a self-consistent manner. 



5 DISCUSSION 

The very small amplitude of the ellipsoidal variations (~ 7.2 • 10"') 
for this period, can only be explained by a body of planetary mass. 
Increasing the mass of the companion would significantly affect 
the amplitude of the ellipsoidal variations (Fig. [3]- bottom), there- 
fore we are confident that the mass of the companion KOI-13.01 is 
8.3±1.25 Mj (i.e a Icr upper limit of 9.55Mj), assuming that there 
are no other contributions to the Kepler light curve that have so far 
been ignored in our analysis. It is worth noting that the ellipsoidal 
component of the light curve is visibly non-symmetric about phase 
0.5 due to Doppler Beaming (Figure[3]l. The difference in heights is 
~ 4.3 ■ 10"^, i.e. 3 per cent of the total phase function signal (assum- 
ing a simple model which neglects reddening). The mass we derive 
for the unseen companion to K0I-13A, places the object below the 
minimum mass for deuterium burning (Spiegel et al. 201 1). 

One possible source of error in our estimation of the compan- 
ion's mass, is the uncertainty in the contribution of the companion 
A-star to the total Kepler light curve. If we increase the contribu- 
tion from the nearby stellar component, and the flux we attribute 
to K0I-13A itself decreases, the mass of the planet becomes cor- 
respondingly larger, as the relative contribution of the ellipsoidal 
variations goes up. However, the effect is rather small when com- 
pared to our other sources of error. Specifically, a factor of ±0.2 
magnitudes, leads to a change in the mass of ±0.05 Afj for the com- 
panion. If the Howell et al. (201 1) delta-magnitude is used, then the 
mass of KOI- 13b will decrease by ~ 20 per cent; this implies the 
body is even more planet-like. 

Another possible source of error in our treatment of the system 
is our assumption that the companion behaves as a Lambert sphere. 
If the correct phase function should be represented by a geomet- 
rical sphere, then the best-fit mass of the system will actually be 
significantly reduced by a systematic factor of 36 per cent, down to 
Mc = 5.3Mj. See IVIislis et al. (2011) for a more thorough discus- 
sion of this issue. For this paper, we are using the Lambert sphere 
phase function, which gives the maximum Mc ■ 

The only real way to increase the mass of the unseen compan- 
ion into the brown-dwarf or stellar regime, is to add yet another un- 
resolved (in the Szabo et al. 201 1 lucky imaging), and comparably 
bright, companion to the Kepler light curve. This seems unlikely, 
given that the spectral analysis of Szabo et al. (2011) is consis- 
tent with the photometry and the model isochrones. Ultimately, the 
mass of KOI-13.01 can only be fully determined with high spa- 
tial resolution (to resolve the components A and B), time-resolved 
spectroscopy. We estimate that the radial-velocity amplitude for 
the KOI-13A system will be A" ~ 870 m s"' for our best-fit mass. 
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Figure 5. Tgq vs /?c for planets found on exoplanet.org. KOI-13.01 is in 
the middle right part of the plot. The solid line is the Rq oc T], -^ model 
from Laughlin et al. (2011). For the plot we have used exoplanets.org 
|http://exoplanets.org/j . 

Although measuring precise radial-velocities for a rapidly rotating 
and spatially-blended A-star is not simple, the formula in Battaglia 
et al. (2008) suggest that it should be eminently achievable, even 
with a small telescope, given sufficient signal-to-noise. Ignoring 
the complications arising from disentangling the two spectra, we 
find that a signal-to-noise ratio of 10, 000 per resolution element, 
should enable us to reach a velocity precision of around 100 me- 
tre/second. 

The low a mplitude of th e ellipsoidal variations rules out the 
suggestion from lSzabcl ( 1201 ll) that the companion could be an M 
star (Fig.IDl. Therefore KOI-13.01 is an exoplanet, and it is now il- 
lustrative to consider the object in the light of other exoplanets, es- 
pecially perhaps WASP-33b which also orbits an A-s tar. In the Tea 
vs /?c diagram (Fig.|5]l, we see that the current radius jBames et al.l 
I20T l|) is in good agreement with the T^q- Laughlin et al. (2011) 
examine the radius anomaly of exoplanets around their best fit 
R oc rj, ( with 7 = 1.4) relation, and we note that KOI-13.01, 
as a candidate for one the hottest planets known, may provide a 
useful constraint on our understanding of planetary structures and 
atmospheres. 

Mazeh et al. (2011) find a mass of 4±2 and 6±3 Mj based 
on ellipsoidal variations and Doppler beaming effect using Rq, R^ 
from Szabo et al. 2011. These values have since been updated by 
Barnes et al. (2011). In our analysis we found that the compan- 
ion is significantly heavier, smaller and closer to the star, than the 
Mazeh et al. 2011 measurement. Also, Shporer et al. (2011) found 
Mq sin(0 = 9.2± 1 . 1 Mj based solely on the Doppler beaming effect 
(also using the Szabo et al. 201 1 mass). The Doppler beaming sig- 
nature in the light curve is roughly an order of magnitude smaller 
than the ellipsoidal variation, as Shporer et al. (2011) show in their 
paper. Assuming an inclination of i=85.9, the planetary mass from 
Shporer et al. is 8.1 < Mc < 10.3. This mass value is significantly 
heavier than the Mazeh et al. value (both teams are using the same 
algorithm and the same system parameters), but agrees with our 
mass value. 



6 CONCLUSIONS 

We present modeling of the KOI-13.01 Kepler light curve, confirm- 
ing the planetary nature of the companion. Our analysis illustrates 



the wealth of information that can be extracted from high precision 
light curves, in the absence of spectroscopy. We find that the out- 
of-eclipse light curve of KOI-13.01 is dominated by two effects. 
The light contributed by the planet is at the same amplitude as the 
ellipsoidal variations induced in the star. 

We are unable to solve for the relative contributions of ther- 
mal and reflected emission from the planet, however we expect the 
planet to be one of the hottest known, given its close proximity 
to an A-star, and its large radius. The equilibrium temperature for 
any zero-albedo companion to KOI-13.01 is = T^jf ^R,l2a = 
2864K, which is rather hotter than calculated for WASP-33b and 
WASP- 12b. We suggest that infrared observations would help to 
disentangle the thermal and reflected components of the light curve. 

By modelling the light curve, we find that the mass of the 
planet is Mc = 8.3 ± 1.25Mj, robust against any assumptions about 
the albedo of the planet. We also find that the planet is orbiting in a 
circular orbit (e = 0.0±0.05). Our results suggest that the density of 
the planet is ~3 times larger than Jupiter's density. This value is not 
surprising. There are much more dense exoplanets, such as X0-3b 
(density ~6.5 times more dense than Jupiter - exoplanet.org) 

We have studied the case that the system is an M-dwarf or 
brown-dwarf eclipsing companion, but neither of these solutions 
can explain the light curve we observe. The expected radial veloc- 
ity amplitude for KOI- 13 system is ^ ~ 870 m s"', which will 
be relatively easy to measure despite the nature of the primary (a 
rapidly rotating A star) and the contamination from the close A-star 
companion. It would be an interesting exercise to search all of the 
Kepler light curves to look for the signature of ellipsoidal variations 
as a planet detection method, even in the absence of transits. 
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